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O ■ Abstract 
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Q«^ . We investigate the two dimensional lattice Gross-Neveu model in large flavor 

C^ . number limit using the domain-wall fermion formulation, as a toy model of 

Mh, lattice QCD. We study nonperturbative behavior of the restoration of chiral 



symmetry of the domain-wall fermions as the extent of the extra dimension 



^ . i^s) is increased to infinity. We find the the parity broken phase (Aoki phase) 

for finite Ng , and study the phase diagram, which is related to the mechanism 

of the chiral restoration in A^^ ^ oo limit. The continuum limit is taken 

and 0{a) scaling violation of observables vanishes in A^^ ^ oo limit. We also 

examine the systematic dependencies of observables to the parameters. 
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I. INTRODUCTION 

Chiral symmetry is one of the important properties to understand the hadron physics 
and the phase transition in the thermodynamics of field theories. Pions are regarded as 
the pseudo Nambu-Goldstone bosons associated with the spontaneous breakdown of chiral 
symmetry. Physics of the phase transition between the confining phase (Hadron phase) and 
the deconfining phase (Quark-Gluon Plasma phase) in QCD is interested from the theoretical 
and experimental point of view. The lattice field theory is one of the most powerful tool for 
such important physics beyond perturbations. 

However there is a problem to define chiral symmetry on a lattice. The definition of 
chiral symmetry on the lattice is one of long-standing problems of the lattice field theory. 
This problem is called as "no-go theorem" ||l]0: the unwanted species of fermions appear 
in chiral symmetric theory on the lattice. To avoid this, chiral symmetry has to be broken 
by adding the Wilson term to the Lagrangian 0,^. This formulation is known as Wilson 
fermion (WF). In order to obtain the chiral symmetric theory in the continuum limit, one 
has to fine tune the quark mass parameter to cancel the additive quantum correction, which 
is a non trivial task in numerical simulations. Besides this fine tuning difficulty, the physical 
prediction from WF has 0{a) scaling violation due to the absence of chiral symmetry, where 
a is the lattice spacing. One must calculate on a fine lattice to get precise predictions in 
WF. 

Few years ago, the domain-wall fermion (DWF) is proposed [^,^ as the new formulation 
of lattice fermion. DWF is formulated as D + 1 dimensional Wilson fermion with the free 
boundary condition in the extra dimension |0 . When the extra dimension becomes infinite, 
DWF has a massless fermion with chiral symmetry. Because of this property, the current 
quark mass term in DWF receives multiplicative renormalization in contrast with WF, and 
this fermion formulation does not need the fine tuning in order to restore chiral symmetry 
on the lattice 0,^. The 0{a) scaling violation is expected to vanish in large Ng limit, which 
means smaller scaling violations than in WF. Thus DWF has desirable property for defining 



the lattice fermion especially for cliiral symmetric or near massless cases. 

The numerical simulations of the domain- wall QCD (DWQCD) have been already carried 
out 0|1^ . As a trade off of the above ideal properties, one needs a large amounts of CPU 
time for computer simulations of DWQCD. Furthermore the new model has larger number 
of parameters than that of WF, and results of simulations have complicated dependencies 
to parameters. 

To clarify the nonperturbative properties of DWF, we examine the formalism for the 
solvable theory, Gross-Neveu model (GN model) in two dimensions, which shares common 
natures with QCD. Our main purpose in this paper is to understand the restoration of 
chiral symmetry. For finite Ng, for which the lattice simulations are performed, we will see 
that chiral symmetry is broken, and the model can be seen as "improved Wilson fermion" . 
DWF for finite A^^ has the parity broken phase similar to WF [jll]-^. We also study the 
continuum limit of DWF model. We will see the 0{a) becomes small as Ns is increased and 
vanishes in the limit A^s ^ oo. 

This paper is organized as follows: In Sec.y, after a brief review of the GN model in the 
continuum, the lattice GN model with DWF is formulated and the effective potentials in 
large flavor (A^) limit are calculate when the extent of extra dimension (A^^) is both infinite 



and finite. In Sec.|T| we calculate the continuum limit in infinite and finite A^s case and 
discuss the restoration of chiral symmetry. We will see that the model has chiral symmetry 
on the lattice for infinite A^^ , while the symmetry is restored only in the continuum limit with 



fine tuning for finite A^^. In Sec.|^ we analyze the structure of the chiral phase boundary 
between the parity symmetric and the parity broken phase (Aoki phase) for the finite lattice 
spacing by solving the gap equation and discuss the necessity of fine tuning to restore the 
chiral symmetry. In Sec.[V] we study the parameter dependences of the lattice observables 
and in Sec.[V| we discuss the way of taking the continuum limit. We will see that the correct 
continuum limit are taken from the lattice model and 0{a) scaling violation vanishes for 
A^s — > oo limit. We conclude with a summary in Sec. |Vli| . 



II. ACTION AND THE EFFECTIVE POTENTIAL 



A. Continuum Gross Neveu model 



The investigation of the GN model is a good test |Tl|JT^-0 for the nonperturbative 



behavior of QCD since the two theories share common properties: the feature of asymptotic 
freedom, chiral symmetry and its spontaneous breakdown. 

The two dimensional continuum Gross-Neveu model in Euclidean space is defined by the 
action 

2 

S = I dx'i^Pij.d, + m)^ - ^m^Y + (^^5^)']} , (1) 

where ip is an A^-component fermion field. The effective potential in large A^ limit in the 
continuum theory is given as 

v;„ = -™. + i-(.^ + tf)to^. (2) 

where m is the renormalized mass and A is the scale parameter. If the momentum integration 
is regularized by a cutoff, M, the renormalization for bare mass mo and the bare coupling 
constant g'^ are 

m = ^ , (3) 

? = 2;^^"a^' (^) 

where the latter shows the asymptotically free. (We set the renormalized coupling constant 
to unity) 

The auxiliary field a and 11 relate to the fermion condensations by equation of motion 

9^ - g"^ - 

a = mo- — V'V', n = —'^Ji-i^ilJ . (5) 

When rriQ vanishes the model shows chiral symmetry, which is expressed by the 0(2) rota- 
tional invariance of effective action (^ in (a, 11) space. The stationary point of the effective 
action (^ satisfies cr^ + 11^ = A^, which manifests the spontaneous breakdown of chiral 
symmetry. 



B. Lattice model with DWF 

DWF is formulated as Wilson fermion (WF) in D + 1 dimensions, or equivalently Ns 
flavor WF with a flavor mixing, which has negative Wilson term obeying the free boundary 
condition on the edges in the extra dimension 0. 

The action of DWF is given as follows: 



^free = «' E E ^(^^ s)Df'''^'{m, S; 71, t)^(n, t) , (6) 

s,t m,n 

where 

2 

D^''^''{m, s; n, t) = E (yiiC^,{m, n)5s,t - W{m, n)Ss,t + PR^s+u^m^n + Ph^s^t+i^m^n, (7) 

C^{m, n) = -^ [5rn+ii,n - 5rn-il,n] , (8) 

Td 

W{m,n) = {I- M)6m,n + -^ [2(5„^,n - S^+f,,n " Sm,n+il] ■ (9) 

Pr/l = (1 ±75)72 are the projection operators into the right- and the left-handed mode, 
a and a^ are the lattice spacing in two dimensions and the third dimensions respectively. 
(T^'s are defined as ai = i and o"2 = 1 in two dimensions, s and t are the indecies of extra 
dimension with 1 < s,t < A^^. Here, r(> 0) is Wilson coupling constant and M is the 
domain- wall mass height (DW-mass). The boundary condition in the third direction takes 
the free boundary condition. In the foUowings, we take a = ag and r = 1 for simplicity. 

If one sets < M < 2, there is single light Dirac fermion in the spectrum of this free 
action whose right (left)-handed components stays near s = l{Ns), 

q{n) = Pn^Pin, s = I) + Pl^u, s = N,) . (10) 

The mass of this light quark, g(n), is exponentially suppressed for large Ng, 

rriga ~ (1 - M)^^ . (11) 

Wilson term in the action avoids the species doubling problem. The doublers and other 
A^s — 1 bulk fermions acquire the cut-off order mass and decouple from low energy physics. 



As a toy model of the lattice DWQCD, we define the two dimensional lattice GN model 
with A^ flavors: 



S = Sfree + o^ Yl cimfq{n)q{n) - a^ ^ 



^ {q{n)q{n)Y + ^ {q{,n)i-i^q{n)Y 



:i2) 



where we abbreviate ■ip{n,s)tp{n,t) = J2iLi'4' {n,s)ijj\n,t). ruf is the "current quark mass" 
in order to give the mass to the fermion[|. In perturbation of DWQCD, rrif receives the 
multiplicative renormalization in A^^ ~^ oo case 0,§|. DW-mass M, on the other hand, 
receives the additive renormalization 0,§], because M corresponds to Wilson mass term. 
We will see that two different couplings, g"^ and g"^, are needed for chiral symmetry in 
general. 

The action (|12|) can be rewritten into the following equivalent action using the auxiliary 
fields crln) and Il(n): 



S = Sfree + ^^ ^ q{n) {aa{n) + i'y^aUin)} q{n) 



+ «^E 



A^ 



A^ 



77-^ Win) -rUf} +—^U{n 



2gy ' ' ^' 2gl- 

The auxiliary fields are related to condensations of fermion; 



(13) 



a{n) = mf- -^q{n)q{n) , U{n) 



9l 

N 



q{n)i-i^q{n) 



(14) 



from the equations of motion. 



C. the effective potential 

We calculate the effective potential of the two dimensional lattice GN model in large A^ 
limit. In large A^ limit, in which the quantum fluctuation of <y{n) and 11 (ra) is suppressed 
and the mean field approximation, (j{n) -^ a and J\{n) -^ 11, becomes exact, the effective 



^ The parameters, M, r and mj, in this paper are opposite sign with the Ref. [18 



potential of GN model is obtained by exponentiating the fermion determinant which is 
calculated by the integration of the fermion fields ip: 



,.^^^^ -a^Et^C— /)^+4n^l 



detD(a,n)e ^^ ^ ^ ^ = ^-a v^ff 



(15) 



This determinant can be calculated by employing the technique of the propagator matrix 
1^,20|. The action introduced in the previous subsection is rewritten by the matrix repre- 



sentation, in the momentum space. From ([TsD, 



s,t m,n 



C^ -w 



\ 



+ 



(oA 



v»V 



5s,t+i + 



( , o^ 



auj 



^s^N^tA + 



( Q acut^ 



-W -C 
( 



fnA 







SsaS, 



t,Ns 



J 



Ss.t + 



\ 







ylOj 



K 



s+l 



i/jL{n,t) 



(16) 



where C^ = sin(pia) + i sin(p2fl), uj = a + ill and W = {1 — M) + I]^=i[l ~ cos(p^a)]. After 
some calculations similar to Ref. [|I^,^], we have 



detD(a,n) =det 



' acu 0^ 



1 



T 



-Ns 



^ 0^ 



V 



auo^ 



(17) 



where T is the transfer matrix along the extra dimension defined by 



T = e 



asHs 



( 1 
W 



1 



c 



\ 



V 



-c'^ , w + c^^c J 



J.-1 



( c^ic + w , c^\ 



V 



w 



w 
l_ 

W J 



:i8) 



Hs is the Hamiltonian along to the extra dimension. The overall factor in Eq. (^) is 
omitted. Diagonalized T~\ whose eigenvalues are A and 1/A, we obtain a explicit formula 
for the determinant, 



det D{a, n) = n [f{M, iV„ p^)a\a^ + H^) + G(M, iV„ p^)aa + H{M, N^, p^) 
where 



(19) 



Pi^ 



H{M,N,,p,) = — [-(A^= - A-^0(2 -h) + (A^» + A-^^/ 



Ns 



-N, 



\-''^){2-h) + {X''^ + X-''')f 



1 



A 



W 



2W 



[h + f] (|A|>1 



1 1 



A 2iy 



[l-M) + J2[l-cos{p,a)] 



p 



[h-f] (|A|<1), 

sin^(pia) + sin^(p2a-) 



(20) 

(21) 

(22) 
(23) 
(24) 

(25) 



Substituting (19) into (15), we obtain the effective potential of the two dimensional lattice 
GN model in the large N limit: 

1 , .9 1 



K 



eff 



I{a,U,M,N, 



_1_ (^ _ raff + ^U' - I(a, U, M, N,) 

n/a (fp 



29V^ -^" ' 2gr 



In 



Fa^(a^ + W) + Ga(j + H 



(26) 

(27) 



r/a (27r)2 

This effective potential is symmetric under M ^ 6 — M. 

Here we comment about the absence of the bosonic fields (Pauli-Villars fields) in our 
model which is employed in (full) QCD [|7|,|9|,p!0|l . Since the gauge field in QCD equally 
interacts with not only light quark field, q{n), but also heavy (bulk) fermions ip{n,s), one 
should introduce the Pauli-Villars boson to subtract this heavy fermion effects. In GN 
model a and 11 play analogous role to the gauge field in full QCD except that they couples 
with q{n) only. Thus it is not necessary to introduce the subtraction in our model. One 
the other hand one could think about a new four Fermi interaction model, in which whole 
ipi^ri, s),s = 1,- ■ ■ ,Ns equally couple to auxiliary fields, a and 11. However these fields are 
constituted from both light and heavy fermions and don't show the chiral property in such 
a new model. 



III. CHIRAL SYMMETRY RESTORATION 

First we show how chiral symmetry restores in this model in the case of both Ng = oo 
and A^s < oo by examining the effective potential given in previous section. 



8 



We will see that for infinite Ns case chiral symmetry is exact even for finite lattice 
spacing, a > 0, without fine tuning for bare mass parameters. The situation for finite A^^ 
case, on the other hand, is much like that of Wilson fermion action. The continuum limit 
has to be taken, and at same time, the bare mass parameter must be tuned finely for chiral 
symmetric effective potential for Ng < oo. 

A. The effective potential in 'W^ = oo" case 

In this subsection we calculate the expression of the effective potential for Ng = oo case. 

For large A^^, one easily sees that the dominant contributions for function I {a, U) in the 
effective potential (^) are the functions "F" and "if", which behave as A^" in A^^ — i> cx3 
limit, from Eqs. ( pO|) and (^). The chiral breaking term "Gaa" in (|27|) can be ignore in 



the limit Ng -^ oo and (^) could be written as 



J(a, n) = r^^ 4^Mh + Fa\a' + H^)] . (28) 

So the effective potential (p6|) is a function of a"^ + IP, which is invariant under the 0(2) 
rotation if g"^ = g^. We emphasize that this is chiral symmetry even before taking the 
continuum limit, a ^ 0. 

The continuum limit of the effective potential can be evaluated by separating divergent 
parts and finite parts from (EBp. Rewriting (P5|) into a integration form 



/(^,n)=/ dpK{p) , K{p)= ^^,-rr , (29) 

one can pick up the the divergent part in a — i> 0, near zero fermion momentum, p^a = (0, 0), 
for pa"^ ~ 0. Since the divergent part of K{p) behaves as 

H I 



Fa? M^{2-Myp^ 
in a — i> limit, the function K{p) becomes 



EpJ. (30) 



^'^>-/t(g^ /^gl^-.P + ^°'^-^>- '^'' 



9 



where fu = M{2- M) and 

C„(M,A.,)=r^Z|5f4zl, (32) 

with ^^ = p^a. The factor /m appears as the normahzation factor of the propagator: 
(qq) ~ fMiiPtilti)"^- The wavefunction of the massless eigenmode has finite width in the s 
direction, and /m is the ratio of q{n) to the the zero mode. 

In A^s = oo case, when < M < 2, only the momentum around p^a = (0, 0) dominates 
in ([28| ) and the contributions of doublers, p^a = (vr, 0), (0, vr) and (vr,7r), are removed com- 
pletely. This means that the doublers decouple from the physical spectrum. For 2 < M < 4 
the momenta p^a = (vr, 0) and (0, vr) become physical poles in the momentum integral, while 
the remaining mode at p^a = (vr, vr) is dominant for 4 < M < 6. In these two regions of M, 
the normalization factor Jm becomes (M — 2) (4 — M) and (M — 4) (6 — M) respectively. For 
M < and M > 6, the no physical pole emerges. 

By evaluation of the first term in the right hand side in Eq. (|3l|) , we find 

K{p) = ^ In -^ + Co{M, Ns), (33) 

47r a^JMP 

where the new constant Cq is defined by Cq = Cq + C'q with 
By substituting this expression into Eq. (^), we obtain 

/ = JA (^2 ^ j^2\ 1^ O^^fll (^' + n^) ^ (j^ / 2 ^ n2) . (35) 

4-/]- V / e ^ ^ 

Therefore the continuum limit of the effective potential for N^ = oo case is 

We find that the two coupling constants can be same as each other g^ = g1 = gl from above 
result, which is in contrast to the GN model with WF. As mentioned above, since the effective 
potential has chiral symmetry for finite lattice spacing, the fine tuning is unnecessary for 
chiral symmetric continuum limit. 

10 



After redefining for a and 11 such as a^ = JmCi and TIr = /mH, we renormalize the 
coupling constant g and m/ as follows: 



m 



^f „_ / 2 ^ ^ JM ,._ 1 



where A is the scale parameter and m is the renormalized mass parameter. One realizes that 
the current quark mass term, m/, receives the multiplicative renormalization (^), similar 
to the perturbation theory in infinitely large Ng. With this choice of scaling relations the 
continuum limit of the effective potential for A^, = cxd is finally given as 

K// = -man + ^ (4 + ^l) In ^^^, (39) 

which is identical with the continuum theory (Q). 

B. Chiral restoration in the continuum limit (A^ finite case) 

As we have seen in previous subsection, the limit Ng -^ oo guarantees chiral symmetry 
on a lattice. Nevertheless, to understand the behavior of finite Ns theory is important for 
lattice QCD simulations using DWF. 

Let us start the finite A^, analysis with the effective potential (^^- The chiral breaking 
term, Gaa, in (^) has the important role for finite A^, and the GN model does not have 
chiral symmetry any more. Chiral symmetry can be restored in the continuum limit with fine 
tuning. In this sense, finite A^s model is similar to WF. In order to restore the chiral symmetry 
without fine tuning, one should take the limit Ns —>■ oo before taking the continuum limit 
to obtain chiral symmetry on lattice. 

The continuum limit for finite A^^ can be taken following the same procedures in the 
previous subsection. The difference from the calculation in A^^ — >^ oo is a necessity to shift 
the auxiliary field: 

a' = cx- ^ '— . (40 

a 

11 



The effective potential becomes as follows: 



K 



eff 



'^gl 



a 



ruf 



;i-M) 



Ns 



'^91 



tf- j(a,n) 



j(a,n) 



ir/a (fp 



n/a (27r)2 



In 



H' + Fa'{a" + W) + G'aa' 



(41) 
(42) 



where 



Ns 



H' = H + F{1 - My^^ + G{1 - My^ , G' = G + 2F(1 - M) 



Expanding /(o", 11) into a power series of lattice spacing: 



\Ns 



(43) 



/ = Jo + /l + /2 + 






In 



Va (27r)2 

(-1)" /•'^/^ dPp 



H' + Fa\a''' + 11^) 
G'aa' 



n J~Tv/a (27r)^ 



(n> 1) 



_/7' + Fa2(a'2 + n2)_ 

One can see that Jq preserve 0(2) rotation in (o"', 11) plane, while /„(n > 1) break it. 
Jo is calculated in the same way as in iV^ = oo case: 



(44) 
(45) 

(46) 



h = / dpK{p) , K{p) = / 



la (27r)2 H^ 



(47) 



Fa 



7 + P 



Here we take the divergent part in a — *> 0. The shifted function FL' has similar behavior to 



the limit Ng ^^ oo: 



H' 



with 



/ 



M 



JM 7 .Pfii 



M(2 - M) 



(48) 



1 - (1 - M)2^= • 
Therefore the logarithmic divergent term has similar form as continuum theory: 



(49) 



47r ^ 



'Pm {<y" + n^ 



Co(^'' + tf 



(50) 



where the new constant Cq is defined by Cq = Co + CL with 
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Co(M,iV, 



'^ <fi fu J2f,^fM- -p- 



TT (27r) 



H' 



JM 1^/1 4/j 



(51) 



and 



TT/a d^p 



fM.l^. 



+ CUM,Ns 



(52) 



The coefficient Co as a function of M for various Ns is plotted in Fig.|l[ 

Remaining J„ in the continuum hmit is easily calculated. Ji {I2) is a linear divergent 
(constant) term in a — *> limit, while /„ (n > 3) vanishes: 

a' n d?i G' 



T =-r 

' a ' aJ-.{2'KYH' 



-a'^C2 



-a 



,2i r d^i 






n2 



(53) 
(54) 



2 7-. (27r)2 

Fig-H shows these coefficients as a function of A^^ for various M . 

From above calculations, the effective potential in finite A^^ case in the continuum limit 
is obtained as follows: 



K 



eff 



TTlr 

.9a 



a J f. 



M'"'^wr^'''^''V7r,'"' 



+ 



M~ VIE ''^^'' ""^ e 



where 



rrif = nif 



[1 - M)^= 



0"iJ = /AfO"' 



n 



R 



/AfH. 



(55) 



(56) 



If the coupling constants and the mass parameter are renormalized as 



1 



/: 



M 



1 



2^"^° + ^^-47^^a2A2^ 



1 (m'CA 
jM \9a a J 

we obtain the correct effective potential of the continuum theory 



(57) 
(58) 
(59) 
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We emphasize that the current quark mass term, ruf, receives the 0(1/ a) additive renor- 
mahzation for A^^ = finite case. As we describe previously, fine tuning for bare mass pa- 
rameter in (|59D is needed for finite A^^, and chiral symmetry is restored in the continuum 
hmit. 

The restoration of chiral symmetry in A^^ ^ oo could also be seen in the scaling relation 



( p7H59|) . The coefficients Ci and C2 represent the magnitudes of the explicit breaking of 
chiral symmetry. Ci is the additive mass counter term, while C2 is the miss match between 
the quadratic terms of scalar and pseudo scalar particle by the quantum correction. If 
we restrict DW-mass to < M < 2, as shown in Fig.^ the coefficients Ci and C2 are 
decreased rapidly as Ns becomes large. The effects of shifting of the fields and the additive 
renormalization of the mass parameter in ([56| ) and (|59D rapidly vanishes with increasing Ng, 
thus the necessity of fine tuning becomes absent in A^^ ^ 00 limit. 

On the other hand, even when M is set out of the region (0,2), the scaling relations 
( p7| - |59| ) lead the chiral symmetric continuum limit. In this case Ci and C2 don't vanish in 



large A^s limit and the fine tuning is necessary, similar to WF. 

IV. THE PARITY BROKEN PHASE AND THE RESTORATION OF CHIRAL 

SYMMETRY 

For lattice QCD with WF, the existence of massless pion for finite lattice spacing is ex- 



plained by the parity broken phase picture proposed by Aoki UTTj]. Although chiral symmetry 
is explicitly broken in WF, one can tune the mass parameter to obtain an exact massless 
pion in the spectrum even for finite lattice spacing. This cannot be understood by the ordi- 
nary picture of Nambu-Goldstone boson in the continuum theory, in which chiral symmetry 
is the exact symmetry of the action and is broken spontaneously. 

Aoki examined the GN model and lattice QCD with WF for a finite lattice spacing, 
and found that the parity symmetric phase and parity (-fiavor) spontaneously broken phase 
coexist in the parameter space of the model. The parity broken phase is characterized by the 
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non-zero condensation of the pseudo scalar density, (11) = {ipi'j^i') 7^ 0. Provided a second- 
order phase transition separating the two phases from each other, pion becomes massless at 
the phase transition point, which is regarded as a massless particle accompanying with the 
continuous phase transition. 

Before starting analysis of phase diagram of DWF model for general Ng, we note the 
equivalence between A'^s = 1 DWF model and Wilson fermion formalism. The effective 
potential of WF, 

1 , .. 1 



V\ 



w 



2,1^^-^^^ ^2,1 



U' 



Va (Pk 



In 



E 



sin^ ki,a 



+ 



("'^IV' 



cos kfj,a) 



+ ffi 



(60) 



can be seen easily by substituting Ng = l,ac = a — {1 — M) in (^61) . The phase boundary of 
Aoki phase forms the three cusps which reach the weak-coupling limit 5f^ = 0atM=l,3,5 
in WF [|lT]-[l^ and A^^ = 1 DWF model. Three cusps arise from the fact that the doublers 
at the conventional continuum limit, {g'^,M) = (0,1), become physical massless modes at 
at M = 3,5. 

On the other hand, in A^^ = 00 limit, the parity broken phase does not exist because of 
the restoration of chiral symmetry, whose explicit breaking causes the parity broken phase. 

For 1 < A^s < 00, it is expected that the parity broken phase exists and the chiral phase 
boundary forms the cusps as surmised from the Wilson like behavior in the previous section. 
The region of the broken phase is distorted and shrinks rapidly as Ng is increased, and 
vanishes in A^^ ~^ 00 limit. 

The parity symmetry is spontaneously broken at parameter point {rrif, g"^; M, Ng), where 
the gap equations have a stable solution 11 7^ 0. Setting the parameters of this model 
{mf^gl^gl^M^Ng) , observables are calculated by solving gap equations, 

dVeff 



da 

dVeff 
dli 



- (a - rrif) - aT{a, H) - -^(a, H) = 



9% 

n 



a 



9l 



^(a,n) 



0, 



(61) 
(62) 
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where 



with ^^ = p^a. In the case of g^ = g1 = g"^, the position of the phase boundary in the 
parameter space could be obtained by solving these gap equations for 11 = e. Taking e — > 
limit, (61) (|6^) define g^ and mf as functions of a which is nothing but the parametric 



representation of the phase boundary if the phase transition is continuous. 

Let us note that DWGN model with g^ = g^ = g^ has only second order phase transition 
between the parity symmetric and the broken phases. This is in contrast with WF case [l7|, 
in which theory a first order phase transition for g"^ ^ g\ is found. 



By differentiating (|62D with respect to 11, one can easily check from that pion mass 
squared 

^1 



Ml(xMl= "^ ' ^ - 



flf dW fl 



9 



(65) 



exactly vanishes at the phase boundary. The factor, 1/ flj, in front of the r.h.s corrects the 
normalization of auxiliary field 11 in the same sense as explained in the previous section. 



(See below (H).) 

Before drawing whole phase diagrams in the parameter space {rrif, g"^), let us first examine 
the positions of the three "tips of cusps" , which are the intersection points between g"^ = 
plane and the critical line. These tips correspond to the continuum limits g^ ^ similar 
to WF. The positions can be obtained by the asymptotic form of the integrals in (|6TD and 
(|6^). The divergent behaviors of the integrals around p^a = (vrrri, vrn), with m,n = or 1, 
are 

^(a,n = 0)^ y / -^ = 2 , (66) 

Q{aU = Q-)^ V r d\ -2(l-M + 20^- 
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where Afj, and S^ are functions oil,a,M and Ng. From these expressions it is easy to see that 
the r.h.s of the gap equations (^) ( |62D have a logarithmic divergence at o" = (1 — M + 2/)^" 
with / = 0,1,2. This fact leads that for each / the phase boundary intercepts with (7^ = 
at a point rrif = (1 — M + 2/)^^ Each of three critical points corresponds to the massless 
particle pole of momentum 

p^(/ = 0) = (0,0) , p^(/=l) = (7r/a,0),(0,7r/a) , p^(/ = 2) = (yr/a, vr/a) . (68) 

The positions of three points move as a function of M, which is shown in Fig.|^. For M in 
the region of (2/,2(/ + 1)), the chiral point for momentum mode Pfj,{l) in (|6^) converges to 
rrif = as Ns increases. Thus massless pion can be automatically obtained at ruf = in 
the limit oi g^ = by increasing A^^. Other two critical points move rapidly to |m/| -^ 00 
in A^s — ^ C)0 limit. 

The actual calculation for the chiral phase boundary is done by solving the gap equations 
numerically, and results are shown in Figs.^, ^ ^ and |^. The former two figures are for A^^ = 
even cases and the latter two show Ng = odd cases. 

The schematic diagram of the phase boundary between the parity symmetric and broken 
phase in (m/, g"^) plane for fixed (M, Ng) is drawn in Fig.^ for A^^ =even case and in Fig.|^ for 
Ng =odd. The phase boundaries in both cases show queer shapes which have three intercept 
points, ruf = {1 - M + 2/)^= with / = 0, 1, 2, on ^f^ = line. 

We find that the phase diagram for A^^ =even is different from that for A^^ =odd. In 
Ng =even case, one can analytically verify from the sign-definiteness of integrands of the 
gap equations that the phase boundary exists only for rUf > region. (For the notation 
in Ref. [0 and numerical simulations, the broken phase appear in rrif < region.) On 



the other hand, the phase boundary intersects with rrif = line in A^^ =odd case. If g"^ is 
decreased from finite value to zero, the phase boundary moves from positive rrif region to 
rrif = (1 — M + 2/)^^ at g"^ = 0, which becomes negative in some region of M for A^^ =odd. 
The chiral phase boundaries for A^^ = odd are shown in Figs. ^ and 0. The parity broken 
phase lies across ruf = plane in Ng =odd case. As similar for A'^s =even, this boundary for 
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odd Ns also converges to ruf = plane for larger A^^- 

Another interesting observation from the phase diagram is finite g^ effect for the pion 
mass. The effective quark mass from the free propagator analysis (mean field analysis) is 

TUqa ~ (Me - Mf^ , (69) 

from which one may think that there exists optimal M(= M^), where the quark (and pion) 
becomes massless even for finite A^^- In current model, for finite A^^ the parity broken phase 
boundary stays apart from ruf = plane for all M ii g"^ > 0. This fact indicates that 
there is no massless pion at rrif = for finite A^^ and finite g"^. We plot the value of rrif 
where M^ = as a function of M in Fig. [1^. (This is the distance of the phase boundary 
point from the rrif = line. This quantity corresponds to the critical value of the inverse 
hopping parameter, 1/kc, in WF.) We find that for fixed finite A^^ and finite g'^, pion does not 
become massless at ruf = . We realize that this discrepancy from the mean field picture 
in ( |69D is due to the part of the additive mass term which is proportional to Ci in (0). For 
strong coupling, for example g^ = 5.0, rriflM^^ = 0) is almost fiat at around M ~ 2.0. This 
indicates a difficulty of detecting "allowed region" of M by observing pion mass for strong 
coupling region. For g"^ = the theory turns out to be free theory, and M = 1 gives M^ = 
at rrif = 0. If A^^ is increased, mflM.,^ = 0) becomes exponentially small convergence to zero 
no matter whether g^ is finite or (See Fig.|ll|). 



The dependence of the phase diagram to A^^ could be seen in Fig.[T^. In the figure 
for M = 0.9, the cusp of smaller rrif corresponds to the chiral continuum limit of the 
conventional momentum mode p^ = (0, 0) and another cusps show that of the doubler 
modes. The first cusp converges to mj = line while the other cusps diverge to ruf — >■ cxo 
in large A^^ limit. Below a critical coupling, g"^, in Figs.|]and ^, a pair of critical lines forms 
either of three continuum limits (three cusps) similar to WF, which can be understood as 
the signal for the recovery of chiral symmetry. For g"^ > g1 the parity broken phases are 
merged to a uniform structure. We find that g1 increases exponentially with A^^; 

gl-e'^''', (c>0). (70) 



At the same time, the phase boundary exponentially approaches to nif = plane with 
increasing Ng. For any g"^, the phase boundary converge at rrif = plane, which can be 
easily seen from (pl| ) and (0). The width of the parity broken phase scales proportional to 
a^ |21|, and shrinks exponentially for A^s -^ oo . these facts are compatible with the exact 
chiral symmetry for finite lattice spacing in A^^ — i> oo limit, as discussed in previous section. 
On the other hand, at M = —0.1 the chiral continuum limit and the phase boundary 
go away from ruf = line with increasing Ns. This behavior shows the violation of chiral 
symmetry at rUf = even in A^s -^ oo limit for M < 0. 

V. (M, Ns) DEPENDENCIES OF LATTICE OBSERVABLES 

Let us turns to discuss about the physical observables in this model. We choose a, (qq) 
and M^ as physical observables. 

In Fig.|l3| we plot M^ and (qq) = rrif — a as a function of rrif for several values of M 
with g"^ = 1, Ns = 4. We find that there exists finite rrif region in which M^^ is zeroQ. This 
region is nothing but the parity broken phase. The ruf dependence of the {qq) shows a 
discontinuous leap in between the parity broken phase. 



The systematic dependences of observables to the parameter M are shown in Figs.^ 



and |T5[ In Fig.|l4| we plot (qq) as a function of M for several rrif with g^ = 1 and A^s = 20. 
(qq) continuously gains in magnitude with increasing M around M ~ 1. (qq) for rrif = 
turns out to have minimum magnitude around M ~ 2. 

M^ also has the systematic dependence to M. The ratio M^/M^ is a smooth function 



of M, A^s and M^ shows essential characteristics of the pion mass. In Fig.|T5| M^fl^ 



^ We set n = throughout our calculation, and the value of physical observables in the parity 
broken phase is not exact. For example, M^ should not be zero in the parity broken phase. Note 



that all figures after this section except Fig. 13 are obtained in parity symmetric parameter region, 
in which our results are precise. 
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(fVeff/diP is plotted against M. Comparing results of strong coupling {g'^ = 5.0) with 
that of weak coupling [g"^ = 1.0) in this figure, we find that the depression of M^flj near 
M ~ 1 for weak coupling is not very manifest for strong coupling especially in small A^^. 
In fact the minimum of M^flj places at M > 2 for A^^ = 4. This implies that in order to 
find the allowed region of M for the chiral continuum limit by observing the depression of 
pion mass, one needs larger Ng for strong coupling than that needed for weak coupling in 
DWQCD simulations. Note that such allowed region of M is unknown a priori in QCD by 
the additive quantum corrections. 

As we pointed out in previous section, pion mass does not vanish for all M at rrif = 
for finite A^s if 9^ > 0. For A^^ = 4, M^ at rrif = has its minimum at around M = 1.4. For 
larger A'^s, M^{mf = 0) in the region < M < 2 tends to be flat with smaller pion mass. 
From this figures we conclude that chiral symmetry is restored for large A^^ if M is set in 
the region < M < 2. There is another region 2 < M < 4, where pion would be massless 
particle in A^^ ^ oo limit. This region corresponds to the region of massless "pion", which 
is made of the chiral modes at p^a = (tt, 0), (0, tt). 

VI. CONTINUUM LIMIT AND DISAPPEARANCE OF 0{a) SCALING 

VIOLATIONS 

Toward the continuum limit a ^ 0, the lattice bare parameters are tuned according to 
the scaling relation (|57| - ^9D for finite A^^. We plot a and M^ as a function of the lattice 



Aa(M, Ns) = exp 



-—{ — -CoiM,N, 



(71) 



spacing aA in Fig. ^for several A^^ with fixed DW-mass, M ~ 1. 

The lattice bare observables systematically depend on M as shown in previous section. 
From (|58D the lattice spacing (or the scale parameter of the theory) is also a function of M, 

2n f 1 

See Fig.^ for (M, A^s) dependence of Cq. The M dependencies of observables cancel with 
that of the lattice spacing, and the correct continuum value are reproduced at a ^ limit. 
(See Fig.|T|.) 
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The interesting observation is the disappearance of 0{a) scahng violation for large Ns 
limit. From theoretical point of view, the exact chiral symmetry in A^^ ~^ oo limit is expected 
to prohibits the dimension {D -\- 1) operators in the quantum correction, which cause 0{a) 
scaling violation. The slope at aA = curve in Fig. |16| is finite at A^s = 2 while the curve 
at Ng = 20 is fiat near aA = 0. This shows that 0{a) scaling violation vanishes in large A^^ 
limit and the scaling violation proportional to a^ exists. 

Fig.0 shows that the remaining 0{a^) error is small for M ~ 1 in large Ns, which is 
less than a few percentage for aA < 0.5 in this model. The reason why this 0{a'^) scaling 
violation is small near M = 1 could be understood by expanding the inverse integrand of 
the function "Jo(o", 11)" for small a: 

^ = pVflf + cv'tfua' + dv'ifua' + ■ • • (72) 

Since I/Za/ is minimum at M = 1 so the 0{a^^ deviation from the continuum formula is also 
minimum for M = 1. This feature might be similar for QCD simulation, except M receives 
the additive renormalization from the quantum fluctuation of the gauge field. 

The renormalization formula employed above for finite Ng is similar to that of WF. The 
bare mass needs to be fine tuned toward a certain critical value. One the other hand, we 
consider to apply the scaling relation for Ng —* oo (without fine tuning), 

_L_c„ + C.^§ln-4„ (73) 

al = al (74) 

f5-". P5) 

to finite Ns lattice observables. This is similar to what is done in the DWQCD simulation 
in a sense. The result of this calculation for finite A^^ is shown in Fig.^. For each A^^, M^ is 
apt to go to the correct continuum value for large aA but tends to diverge for smaller lattice 
spacing. Such divergent behavior is never seen in current DWQCD simulations @,0 since 
aAqcD is larger than 0.1. 

From ( ^9l) the renormalized mass is expressed as a function of the bare quark mass mf: 
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TTia 



1 U^a-(l-M)^^+/Ci(M,iV,)l . (76) 

Jm9 '- -' 

For each g"^ (and aA), if Ng fulfills a condition; 

N,> N'^ s.t. |am/| > I - (1 - Mf^ + g^Ci{M, N^)\ , (77) 

the renormalized mass approximately becomes that in A^s ^ oo case. Such A^^ can exist for 
finite m/ if < M < 2 because both (1 - M)^= and Ci(M, A^^) go to zero for large A^^- This 
leads that if A^^ is larger than A^^, physical predictions from DWF is saturated as function 
of Ns and could be regarded as the value of A^^ -^ oo. A^^ is a function oi g'^,M and rrif, 
and tends to be larger for smaller rrif. Ng dependence of M^ could be seen in Fig.|19|. For 



large magnitude of mj the A^^ dependence is saturated up to A^^ ~ 12, while M^ varies as a 
function of A^^ near rrif = 0. From this figure for {g'^,M) = (1,0.4), we can estimate A^^ is 
near 10 for \mf\ ~ 0.1 and A^^ ~ 20 for \mf\ ~ 0.02. In this model A^^ goes to small value 
for M ~ 1. For M = 0.9 the value of M^ for A^"^ = 12 is nearly identical to that of A^'^ = 20 
for almost all the region of rrif in Fig.lT^. 



VII. CONCLUSIONS AND DISCUSSIONS 

We have investigated the two dimensional lattice GN model with DWF in large flavor(A^) 
limit, as the toy model of lattice QCD with DWF. By calculating the effective potential we 
study the nonperturbative prospects of this model, which are expected to be qualitatively 
similar to the DWQCD. 

In infinite A^^ case, the effective potential has exact chiral symmetry even for finite lattice 
spacing. The chiral phase boundary places exactly on rrif = line, which shows the fine 
tuning of the mass parameter, ?7i/, becomes needless. The parity broken phase does not 
exist for all coupling constants. Thus the model for Ng = oo has similar properties as 
the continuum theory especially for chiral symmetry, by which the massless pion could be 
understand as NG-boson accompanying with the spontaneous breakdown of the symmetry. 
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In finite Ns case, for which numerical simulations are carried out, is practically important. 
Chiral symmetry is explicitly broken by the finite A^^ effect, that causes the parity broken 
phase with {D + 1) cusps near g"^ = similar to Aoki phase of WF. The restoration of chiral 
symmetry only occurs in the continuum limit with fine tuning of rrif to its critical point, 
which is on the phase boundary of the parity broken phase. By increasing A^^ for fixed a, 
the phase boundary exponentially approaches to ruf = plane in the parameter space and 
the parity broken phase vanishes in A^s ^ oo. If one take the limit A^s — > oo prior to a -^ 
limit, chiral symmetry is restored. 

Another interesting observation is the A^^ dependence of the critical coupling, g"^. For 
g^ < g1, the parity broken phase splits into {D + 1) regions, each of which corresponds 
to a chiral symmetric continuum limit. The restoration of the chiral symmetry is also 
characterized by the fact that g^ exponentially goes to large value with increasing A^^. 

We also show (M, Ng) dependencies of lattice bare observables, which could be some 
hints for DWQCD simulations. When A^^ is finite and g^ > 0, the pion mass at rrif = 
never goes to zero for all M. If M is set into the range (0,2), the pion mass at mj = 
exponentially suppressed with increasing A^^. From the results of mj at M^ = (Fig.|TD|) and 
M^ at nif = (Fig.|T^), we discussed that one needs to take larger A^^ for strong coupling 
than that needed for weak coupling in DWQCD simulations, in order to search the allowed 
region of M for the chiral continuum limit by observing the depression of pion mass. By 
observing the A^^ dependences of the observables we discuss about the criterion, A^^ > A^^, 
for which the physical observables can be considered as approximate values for A^^ -^ oo. 
A^^ is a function oi g'^,M and rrif. 

The observables depend on the value of M even for Ns -^ oo. This dependence cancels 
by the renormalization and the correct continuum theory is obtained. The disappearance 
of 0{a) scaling violation for large A^^ in the continuum limit suggests the probability of 
obtaining the reliable physical predictions for smaller lattice spacing than that in WF. 

Since the GN model in large A^ limit neglects the quantum fiuctuations, and omits the 
gauge fields, we cannot insist that the behavior of lattice QCD with DWF is exactly same 
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as the results in this paper. For example, the DW-mass, M, is shifted into M = M + Const. 
due to the back reaction of the gauge fields. However there was many similarities between 



GN model and QCD for Wilson action ITlHTSl], we expect that the results shown in this 



paper will provide instructive and systematic informations about nonperturbative effects of 
lattice QCD with DWF. 
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FIG. 3. The chiral symmetric continuum limit in {M,amf) plane for iVs=even. The critical 
lines approach to nif = with increasing Ns in each of three region, 21 < M <2{l + 1), / = 0, 1, 2. 
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FIG. 9. Same as Fig.^ for Ns =odd. The parity broken phase lies across mj = line. 
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FIG. 16. a/ A and M^ as a function of aA for fixed M. The continuum limit is taken using 
the (Wilson like) scaling relations in (57-|59|). 



scaling of a 
m/A=0.0 Ns=20 



-0.9 



M=1-1 



M=1.05 



Cont. 



~----1^», 



0-95 



^^n 



0.1 



'^■9 



0.2 0.3 0.4 0.5 

aA 



0.15 



0.14 



0.13 



0.12 



0.11 



scaling of M^ 
m/A=-0.2 Ns=20 



/ / 


/ 

/ 
/ 


tJ\=A.A ' 


Cont. - — - — _ 






-,~~~~-- V,/,„ 



0.1 0.2 0.3 0.4 0.5 

aA 



FIG. 17. cj/A and Ml as a function of a A at iVs = 20. 
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FIG. 18. a/ A and M^ as a function of aA at M = 0.9. The DWF scaling relations in (^|-|75|) 
are employed. 
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FIG. 19. Ns dependence of M^ as a function of mj. 
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